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Abstract 

An elegant quaternionic formulation is given for the Lagrangian advection equation for velocity 
vector potential in fluid dynamics. At first we study the topological significance of a restricted 
conserved quantity viz., stream-helicity and later more realistic configuration of open streamlines 
is figured out. Also, using Clebsch parameterisation of the velocity vector potential yet another 
physical significance for the stream-helicity is provided. Finally we give a Nambu-Poisson formalism 
of the Lagrangian advection equation for velocity vector potential. 
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I. INTRODUCTION 



Quite recently [1], a conserved quantity — stream-helicity — has been conjectured for cer- 
tain cases of incompressible fluid flows. Although Arnold and Khesin [2] have documented 
in their book a rather general topological study of quantities of such type, stream-helicity 
has never been explicitly studied and used in the mathematical fluid dynamics. Its study 
becomes even more important at the face of the discovery that it has been shown to remain 
conserved for at least some kinds of fluid flows. Therefore, an upsurge of motivation to study 
stream-helicity and the dynamics and beautiful mathematical aspects of the velocity vector 
potential (that has got itself involved in defining the stream-helicity) is quite natural. 
In this paper, as a rather interesting product of the study, the possible quaternionic struc- 
ture of the Lagrangian advection equation for the velocity vector potential has been un- 
leashed; readers may be aware that quaternions, a concept originally invented by Hamil- 
ton to generalise complex numbers to M. 4 , obey multiplication rules governed by a simple 
non-commutative division algebra and currently play an important part in the theory of 
4-manifolds, the technical foundation of modern inertial guidance systems in the aerospace 
industry, the control of the orientation of tumbling objects in computer animations etc. 
Also, for the aforementioned advection equation we have given a possible Nambu-Poisson 
formulation. Most importantly, the topological meaning for the stream-helicity has been 
extended to unclosed but tangled streamlines. For the sake of clarity and completeness, we 
have briefly but appropriately summarised included the essence of the paper ([1]) herein in 
the section (IV). 

This paper is organised as follows: Section II deals with the advection equation of the 
velocity vector potential. We give the quaternionic formulation of the advection equation 
in Section III. Section IV is devoted on stream-helicity starting right from its definition 
and Section V discusses the stream-helicity's geometrical and topological interpretations for 
open streamlines. We indroduces Clebsch parameterisation for the vector potential in Sec- 
tion VI. Section VII reformulates the advection equation for the velocity vector potential 
using Nambu brackets. 
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II. ADVECTION EQUATION FOR VELOCITY VECTOR POTENTIAL 



Let us begin with the Euler equation 

| + (h.V)« = -VP (1) 

for three-dimensional, inviscid, incompressible fluid being acted upon by irrotational body 
forces. The pressure P used in the equation (1) includes the effect of such forces also. As we 
have assumed the fluid to be incompressible, i.e., the density is constant, we have taken the 
density to be unity without the loss of generality. The mathematical way of acknowledging 
incompressibility yields for the velocity field u: 

V.w = (2) 

paving way for defining the vector potential £ for the velocity field as follows: 

u = Vx£ (3) 

— * — » 

Obviously, £ is not unique, for, a term V7, 7 being a scalar field can always be added to it 
keeping u unchanged. We shall come back to this issue in the end of Section IV. For the 
time being we put relation (3) in the equation (1) to get: 

^(V x i) + (tf.V)(V x = -VP (4) 

Trivial manipulations using vector algebra suggests: 

[V x (u.V)|]i = e ijk dj{uidi^ k ) 
=>• [V x (u.V)f]i = €ij k {djUi){di^ k ) + eijkUidjdiik 
[(«.V)(V x = [Vx (u.V)(] t - e ijk (d j u l )(d l Z k ) (5) 

Putting relation (5) in the equation (4) we get the following Lagrangian advection equation 
for the velocity vector potential: 

di - - 

+ («.V)f = curl" 1 ?/ (6) 

where fj is defined as: 

rji = e ijk T jk - d t P, (7) 



where we interpret Tj k as a contribution of a stress tensor which, of course, is not symmetric. 
The stress tensor T jk that does all the possible deformations viz., tilt, stretch, rotation etc. 
of the velocity vector potential lines has been defined as 

T lk = (d jUl )(d^ k ) (8) 

Definitely curl _1 f/ is unique only upto a gradient of a scalar, say, 9, an issue taken up again 
in the Section IV. So, we have yet another gauge to be fixed. In the next section, we shall 
note the equation (6), is endowed with a beautiful exposition of quaternionic structure a la 
Gibbon et. al. [3-5] 



III. QUATERNIONIC FORMULATION FOR ADVECTION EQUATION OF £ 

Let us define A = curl rf for the sake of simplicity of writing. Of course, as discussed in 

— * 

the Section II, we need to fix a gauge for A and this possibly can be done as in the many 

standard ways of gauge fixing. Thus, we can rewrite the equation (6) and its Lagrangian 
time derivative respectively as 

— * 

TT = * ( 9a ) 

Dt v ; 

W =P =Dt (9b) 

where p has been appropriately defined. Unlike the vorticity vector in the ideal incompress- 
ible fluid, velocity vector potential is not known to be frozen in the fluid. However, it doesn't 
prevent lines of £ from getting tilted, stretched and rotated due to the fluid motion; rather 
the roller-coaster ride of each line of £ is even more complex than the usually studied field 
lines of the vorticity vector; Cu = V x u. Evidently, A is the relevant vector field which 

— # — # 

brings about the deformation in the £ as the fluid evolves. The information about how £ is 

— * 

deformed is obviously kept hided in the evolution of the relative angle and magnitude of £ 

— * 

with respect to A. Keeping this in mind we define following quantities: 

_ £A _ _ I x A _ £.p _ _ £ x p 

ax = W Xx = ~lT' a ' = Ti' X " = TT (10) 
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Also, if the instantaneous angle between £ and A be A , then 

tan0 A = — =r^— = — (11) 
f.A "a 

=>. A 2 = £ 2 a1; sec 2 A (12) 
Using the relations (10) and (12), following equations may be constructed: 

^(f-A) = £ 2 («a sec 2 A + a p ) (13a) 

^(fxA)=£ 2 x P (13b) 

which can then be modified respectively into the following Lagrangian advection relations 
for a x and Xx- 

Da x _ 2 2 



Dt 

Dxx 



Dt 

The equation (14b) may be rewritten for the magnitude of Xx as 



= Xx - a x + a P ( 14a ) 



-2a x xx + X P (14b) 



DXx 

Dt ~" AV '' x " 



= 1q,\X\ + xl (15) 



where, X ' P = (16) 
Xa 

Again, the equations (9a) and (9b) when recasted, with the help of the definitions (10) and 
their evolution equations (14a) and (14b), to get Lagrangian dynamics for the magnitude of 
£, one arrives at: 

^ = «a£ (17a) 



D\ 



= (xi + (17b) 



Dt 2 

This highlights how significantly the dynamics of velocity vector potential is governed by 
the quantities defined in the relation (10). 

^ ^ — * — * 

Now, let £, Xx and £ x xx be the unit vector along £, % A and £ x x a respectively. One can 



form an orthonormal frame x X\,X\) using them. We wish to study the Lagrangian 
dynamics of the frame. In this frame one can apparently write the vector p in terms of the 
independent unit vectors £, £ x xx an d Xa as 

p = a p i + a(ix X x)+bxx (18) 
=>• X? = 6(1 x Xx) - axx (19) 

We have cross-producted the equation (18) with £ to reach at the relation (19). Therefore, 
obviously a = — (xx-Xp) an d b = £.(x\ x Xp)- Now, defining the Darboux angular velocity 
vector as 

n D = xx + —i (20) 

we can arrive (using the expressions (9a), (10), (14b) and (15)) at the following evolution 
equations determining the Lagrangian dynamics of the frame: 

^ = &d x i (21a) 
^ (| x xa) = x (| x xx) (21b) 



Dt 



= d D x X x (21c) 



Before delving into the quaternionic formulation let us see how complex numbers, precursor 
of quaternions, can enrich the beauty of the equations being dealt with. Lets define 

Ca = «a + ^Xa, Cp = ttp + iXp (22) 

After a little bit of manipulation, putting the expressions (14a) and (15) into use, we can 
arrive at 

^ + Ca"C = (23) 

which can be decomposed into the following two equations 

£|C: 12 



Dt 



+ 2c* a |Ca| 2 -2&[CaCJ =0 (24a) 
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_ -(tan A ) + «Atan' j A + a A + — tan0 A ^ = (24b) 

Dt V "a/ «A 

where the 'overhead bar' denotes the operation of complex conjugation. By the way, by 
picking a function ip, such that (\ = 7757, one can linearise the equation (23) to the form: 



ip Dt 

D 2 i) 



C P V = (25) 



Dt 2 

Mood is now set up for introducing the quaternionic structure of the equations for the 
variables a\, x\, a p , Xp and £. We begin by defining following quaternions as the natural 
tetrads 

Qa = [a x , xxi Q P = [a p , x P \ and Q ? = [0, |] (26) 

In what follows let '<£>' and 'asterisk' denote the Grassman product and the quaternion 
conjugate respectively. Then obviously, 

Qa ® = [a x , Xx] ® [0, 3 = [-&•£ <*xt + Xx * |] (27) 
So, one has by mere inspection 

^ - Qa ® = (28) 

and similarly, 

-^-Q p ®Q 5 = (29) 

Differentiating equation (28) and putting equation (29) into use in the midway of the ele- 
mentary calculations, we get the compatibility condition: 

^± ® Q 5 + Qa ® Qa ® Q € - Q P <8> Q^ = (30) 
which due to associativity yields on taking Grassman product with Q^ 1 from the right 

^ + Qa Qa - Q P = (31) 

that is nothing but the quaternionic Riccati equation that can be linearised using a quater- 
nion such that 



Qa = ^ ® Q; 1 (32) 
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where, Q^ 1 = q^q* is the inverse of Q^,. Using this in the equation (30) and the following 
identity for any arbitrary quaternion Q(x,t): 

«£l = -Q-. 8 (33, 
we arrive at the linearised quaternionic Riccati equation: 



D 2 Q 



Df2 -Q„ Q, (34) 

Now, let us set up a Lagrangian advection equation for Q p analogous to the relation (31). 
To do that let us note that as both xx an d Xp are perpendicular to £, one may write 

ix P = (35) 
Xx x X P + ciX\ + c 2 Xp (36) 



Dx P 



Dt 

where c\ and c 2 are two undetermined scalars. In reaching the expression (36) from the 
relation (35), we have made use of the equation (21a). Using the definition (10) explicitly 
in the relation (36), a bit of manipulation yields 



£x 



D (p\ A 



Dt\i) % 



CiXx + c 2 x P (37) 



^ (jj = <x p j + x i + c 2 x P x i + c 3 | (38) 
where c 3 is yet another unknown scalar. Recalling the relations (10) and (21a), the preceding 



equation (38) when dot-producted with £ yields 

Dct p _ _ 

-j^- = a x a p + xx-Xp + c 3 (39) 

With the information contained in the relations (36) and (39) one can get the coveted 
quaternion equation as follows: 



= Qx ® Q P + Q c (40) 



Dt 

In which, using the identity quaternion (I = [1, 0, 0, 0]) with respect to the Grassman prod- 
uct, we have defined the quaternion Q c as 

Q c = ciQ„ + 4Q A + c£l (41) 

c[, c' 2 and c 3 , in turn, are the redefined unknown scalar given as 

c[ = ci - a x , c 2 = c 2 - a p , and c' 3 I = (c 3 - c x a p - c 2 «a)I + Qa ® Q* p + Q P <S> Q* x (42) 

Thus, unfortunately we could find the Lagrangian quaternionic equation for Q p only at the 
cost of introducing three undetermined scalars. 



IV. STREAM-HELICITY 



The 'stream-helicity' (S) is defined as: 

S= I g.ud 3 x (43) 
Jv 

where V is a volume occupied by the fluid. At this point let us ponder over the aforemen- 
tioned non-uniqueness of the vector potential [6]. For smooth discussion's sake, we assume 
for the time being that the volume is simply connected. Suppose £ — > xi + V7, then from 
the definition (43) of stream-helicity we can find the change SS in S to be: 



5S = V^.ud A x = <h ryu.nd x (44) 

JV JdV 

where n is the unit vector perpendicular to the infinitesimal surface element d 2 x and we have 
used the relation (2) and Gauss divergence theorem. The relation (44) amounts to saying 
that the stream-helicity will be gauge invariant in case the surface dV bounding V, is the 
surface made up of streamlines i.e., u.n = on dV. This condition for gauge invariance is 
rather strong because if u.n 7^ on dV then one cannot seek refuge in Coulomb gauge for 
it is too loosely defined inside V with no information about the outside field whatsoever. 
More starkly, it means that different solenoidal vector potentials inside V can correspond 
to Coulomb potentials of fields which have different structures outside V. Now if we relax 

— * 

the condition that V is simply connected, then the line integrals of £ about the 'holes' in 
the possibly multiply connected region have to be specified in order to have gauge-invariant 
stream-helicity within dV on which u.n = 0. Another good thing about this very gauge is 
that it helps in bypassing the issue of non-uniqueness of curl -1 ?/ in the discussion on the 
conservation of S in this section by causing the R.H.S. of the expression (46) to become 
independent of 9 that naturally infects the second term of the R.H.S. of the equation. 
Let us take total derivative of S w.r.t. time to get: 

Tt = / Di^ d x < 45) 



dS_ 
~dt 



Ji{-VP)d z x + J u.icvxl^rOtfx (46) 



where, D/Dt is the material derivative w.r.t. time and it basically is a shorthand for 
d/dt + u.V '. Again, simple vector algebra suggests: 

(V x ^.(curPV) = V.(|* x curP 1 ^ + 77. |* (47) 



With relation (3) in mind, inserting relation (47) in the equation (46), we have the following: 
-2 J C^Pd 3 x + J V.(?x cur\- 1 rj)d 3 x + J iMe^d^x (48) 



dS 
~dt 



where, equation (2) has been used. If all the terms in the R.H.S. of equation (48) vanish 
then one may set 

§ - < 49 > 

and say that stream-helicity is a conserved quantity. 

[1] picks it up right from here to conjecture that in certain kind of fluid flow the stream- 
helicity is conserved. Further, it has been shown therein that not only does this remain 
conserved but also one can lend stream-helicity a topological interpretation for the closed 
linked stream-tubes and knotted stream-tubes; the stream-helicity is proportional to the 
linkage number in the case of linked stream-tubes for Beltrami vector potential and in the 
case of knotted tubes it is connected to the way in which the knot may be constructed from 
a closed tubular unknotted stream-tube. 

Historically speaking, after almost hundred years of Lord Kelvin's realisation that any linkage 
or any knottedness in the vorticity field at any earlier time should remain conserved at 
all later times in an inviscid and barotropic fluid being acted upon by irrotational body 
forces, Moreau[7] and later Moffatt[8] established an invariant known as helicity which is of 
topological character and encompasses Kelvin's insight. But what Lord Kelvin had missed 
was the possible existence of knotted streamtubes in the steady Euler flows — a fact very 
logically speculated by Moffatt[9]. Actually, inspired by this speculation only, a new quantity 
viz., stream-helicity has been introduced for the case of inviscid and incompressible fluid 
being forced by irrotational body forces. 



V. TOPOLOGICAL MEANING OF S FOR UNCLOSED STREAMLINES 

In this paper we show that this analogy may be stretched further to good use. Just as 
Arnold [10] extended the interpretation the linking of the oriented vortex tubes introduced 
by Moffatt to the situations in which vortex lines are not closed curves but are wound about 
each other infinitely often, we too can identify the stream-helicity as the "asymptotic Hopf 
invariant". Probably the physical situations of the tornado, whirlpool etc. spiraling flows 
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that may cover quite a distance in space can entertain such an interpretation for the stream- 
helicity, for, the existence of the streamlines wound many times around each other is a rather 
strong possibility therein. 

Let us try to be more explicit in the language of differential geometry. For the case of 
incompressible flows being dealt with herein, let the velocity field u be the field on a 3-D 
closed oriented simply connected manifold M with a volume form \x and the phase flow 
of u be 0* : M — > M. Let M has the real cohomology of a three-sphere which means 
if-^M) = H 2 (M) = 0. Then, one knows[ll] a set £ containing exactly one oriented path 
having starting point p and end point q for every pair (p, q) of points of M is always possible 
to choose. As we are concerned with the realistic fluids, in what follows it will suffice for 
us to consider M to be 3-D compact Euclidean domain (M C M 3 ) such that u is tangent 
to the boundary <9M which is, by the way, required for the uniqueness of stream-helicity as 
discussed at the end of the section (II). With these in mind, we consider a pair of points 
x,y G M and select two large numbers T x and T y so that the trajectories ^(x) (t G [0,T X ]) 
and 0*(y) (t G [0,7^]) are closed by joining x and 4> Tx (x) (and similarly, y and <fi Ty (y)) using 
a path belonging to S. Let the two closed curves formed be T x and T y respectively. It is 
assumed that T x and T y do not interest. Obviously, r^. : S 1 — > M and T y : S* 1 — > M are 
smooth mappings of two circumferences to M with disjoint images. In the images, T x (t x ) 
and Ty(t y ) are the velocity vectors found by differentiating the curves with respect to the 
corresponding parameters t x (mod T x ) and t y (mod T y ) respectively. One can now define 
linking number for the closed curves T^S* 1 ) and T^S* 1 ) in M by the Gauss theorem (see 
e.g., [2]) as: 

Mr t)- 1 r r ^^ x )j y (t y ),Y X (t x )-v y (t y )) 

ik(T x ,r y )--j o y o llTx{tx) - Ty (t y w x y m 

where (•, •, •) and || • || are the symbols for vector triple product and the norm respectively. 
Now, the asymptotic linking number for 0*(x) and of u is traditionally defined as: 

X,(x,y)= hm (51) 

The limit is known to be existing almost everywhere. This relation (51) basically the Gauss 
linking number of (^(x) and ^(y) i.e., 

Ux v) ~ (-) hm r C V (Ht x ),y(t y ),x(Q-y(t y )) 
^^-WJt^ooI J \\x(t x ) - y(t y W x tv (M) 
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where x(t) = and y(t) = (^(y) for convenience. Subsequently, the average self-linking 

number of u may be proposed as: 

As = / / Xa(x,y)niH2 (53) 
Jm Jm 

One may now invoke helicity theorem due to Arnold[10] for the case of solenoidal velocity 
field u to seek a topological interpretation for the stream-helicity for the fluid flow having 
unclosed streamlines as: 

S = Xu (54) 

This completes our discussion on the possible topological interpretation for the restrictedly 
conserved quantity stream-helicity for the certain classes of flow. 



VI. CLEBSCH PARAMETRISATION OF £ 

In the form notation, £ is a 1-form such that the Clebsch parameterisation of the velocity 
potential £ using three independent scalars a, f3 and 7 can be represented as: 

£ = ad(3 + d-f (55) 

The parameterisation is attainable locally in space owing to the Darboux's theorem. Of 
course, the uniqueness of a and (5 is at stake because the expression (55) keeps its form 
intact after any canonical transformation of a and f3. In the standard vectorial notation, we 
have from (55): 

I* = «V/3 + V7 (56) 
4u = Vx( = VaxV(3 (57) 

The expression (57) regales itself with the standard interpretation that each streamline is 
the line of intersection of a pair of surfaces (rather stream-surfaces) given by a = constant 
and P = constant. Now, the relations (56) and (57) can be dotted to get the following form 
for the 'stream-helicity density' s: 

s = £.m = V.(7m) (58) 

This expression suggests that for the non-vanishing stream-helicity and appropriately well- 
behaved velocity field, 7 must be singular either on the surface at infinity bounding the 
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integration region or in the finite volume of the integration domain. Therefore, the stream- 
helicity becomes an integral over the surface bounding the singularities: 

/ Cud 3 x= / {{-fu).n}d 2 x (59) 

JV JdV 

As a result, a physical interpretation of the stream- helicity is that it is the measure of the 
flux of the velocity field, weighted by 7, through the surfaces that enclose the singularities 
of 7. 

Again from the relations (57) and (58), we can see that 

(u.V)a = 0; (u.V)P = 0; (u.V)7 = s (60) 
which dictates that the Jacobian J of the transformation (assumed invertible) 

x = X(a,0,'y,t) (61) 

is given by: 

j = d(x,y,z) = 1 = 1 ^ 

' d{a,P,j) {(Va x V/3).V 7 } s 

— * 

Thus, the stream-helicity density s is the Jacobian for the inverse transformation of X in 
(61). 



VII. NAMBU STRUCTURE AND HAMILTON! AN FORMULATION 

In 1973, in a remarkable paper, Nambu [12] introduced a generalization of classical Hamil- 
tonian mechanics which is now called Nambu mechanics. In his formalism, he replaced the 
canonical pair of variables found in Hamiltonian mechanics with a triplet of coordinates in 
an odd dimensional space. He formulated his dynamics in terms of ternary operation, the 
Nambu bracket, as opposed to the usual binary Poisson bracket. The original motivation of 
Nambu was to show that Hamiltonian mechanics is not the only formulation that makes a 
statistical mechanics possible. 

A manifold M is called a Nambu-Poisson manifold [13] if there exits a M-multi-linear map 

{,...,} : [C°°(M)f n -> C°°(M). (63) 

This is called Nambu-Poisson bracket of order n V/i, /2, • • • , /Wi-i € C°°(M). This bracket 
satisfies 
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1. {/!,...,/„} = (-l) e W{/o(l),...,/a(n)}, 



2- {/1/2, /3, • • • , /n+l} — /l{/2, f3, ■ ■ ■ , fn+l} + {/l, f3, ■ ■ ■ , fn+l}f: 



2- 



3. It satisfies fundamental identity, 





H 1" {/id • • • j /2n-2, {/l, • • • > /ra-lj /2n-l}} — {/l, • • • 5 /ra-lj {/n, • • 




where a G S n — the symmetric group of n elements — and e(<r) is its parity. Let M denote 
a smooth n-dimensional manifold and C°°(M) the algebra of infinitely differentiable real 
valued functions on M. 

It is known that the Nambu dynamics on a Nambu-Poisson phase space involves n — 1 
so-called Nambu-Hamiltonians H 1: . . . , H n _i e C°°(M) and is governed by the following 
equations of motion 



A solution to the Nambu-Hamilton equations of motion produces an evolution operator 
Ut which by virtue of the fundamental identity preserves the Nambu bracket structure on 
C°°(M). A function / e C°°(M) is a first integral of A jH - li ...^ n _ 1 if and only if 



We focus on trilinear antisymmetric Nambu-Poisson bracket {f,H 1 ,H 2 }, where Hi and 
H2 are a pair of Hamiltonians. Since {/, Hi} represents the contraction of {f, Hi, H2} 
with H 2 , the second Hamiltonian H 2 is always a Casimir of the original Poisson bracket, 
i.e., {/, H 2 } = for any /. It has been found [14, 15] that two and three dimensional 
Euler equations of incompressible fluid can be manifested in terms of Nambu-dynamics. In 
particular, Nambu mechanics is extended to incompressible ideal hydrodynamical flow using 
energy and helicity in 3D. It is equivalent to 




(64) 



{/, Hi, H 2 , ■ ■ ■ , H n _i} = 0. 



dt 



{f,h,H} 



(65) 



where H = | J u-ud 3 x is the Hamiltonian of the fluid with velocity u and vorticity w = Vxm; 
h = \ f u ■ u d 3 x is the helicity. 
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In our case stream-helicity does not seem to participate in the construction of Nambu- 
dynamics by acting like a second Hamiltonian. Structurally the Lagrangian advection 
equation for velocity vector potential (6) without i] term is very much similar to the in- 
compressible 2D hydrodynamics equation. The only difference is that the equation (6) is a 
three dimensional one. But this makes an enormous difference in studying Nambu-Poisson 
formulation of equation (6). 

We overcome this problem by introducing a pair of stream functions. Thus we define 

u = V^i x V^ 2 , (66) 

where ip\ and ip2 are stream functions. One can easily see that V • u — 0. This enables us 
to express the equation (6) as 

^ + {e,^i,^} = curl- 1 r/ ! (67) 

where 

We assume that the derivatives of the two Hamiltonians H\ and if 2 with respect to £ are 



given by = ipi and = 1^2 respectively. We propose that the equation (67) could be 
written as 



<9£ K .( r 5H 1 5H 2 



a +A Tirirr currl? (68) 

where M is a trilinear antisymmetric bracket, which is defined earlier. 



VIII. CONCLUSION 



In the light of the fact that the stream-helicity can turn out to be a well-defined conserved 
quantity in certain kinds of fluid flow, the velocity vector potential, that basically is thought 
to be a rather unphysical quantity, gets a life of its own. This is so because the velocity vector 
potential can interplay with the fluid velocity to construct the restrictedly conserved quantity 
- stream-helicity — which in itself is very much physical as has been demonstrated in this 
paper by highlighting very sound topological meaning for it. Quite possibly, this conservation 
may have kept hidden within itself a symmetry that may be an interesting prospect to 
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study in future. This paper's naive attempt has been to initiate serious study of velocity 
vector potential's dynamics in fluids. It has been showcased that this vector potential can 
regale itself with the mathematical structures galore viz., quaternionic formulation, Clebsch 
parameterisation and Nambu structure. To conclude, herein has been documented probably 
the first coherent, though diverse, study of the dynamics of the velocity vector potential in 
a 3D inviscid incompressible flow; hopefully more works will be done in this direction will 
follow soon. 
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